A REMARKABLE SEQUENCE OF INTEGERS 
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Abstract. A survey of properties of a sequence of coefficients appearing in 
the evaluation of a quartic definite integral is presented. These properties are 
of analytical, combinatorial and number-theoretical nature. 



1. A QUARTIC INTEGRAL 

The problem of explicit evaluation of definite integrals has been greatly simplified 
due to the advances in symbolic languages like Mathematica and Maple. Some years 
ago the first author described in [2^ how he got interested in these topics and the 
appearance of the sequence of rational numbers 

2m — 2k\ / m + k\ fk 



(1.1) d,.„ = 2-^™> 2 , , ,, ,, , 

^ — ' \ m — k / \ m J \ I 

k=i ^ ^ ^ ^ ^ 

for < I < m. These are rational numbers with a simple denominator. The 
numbers 2'^"^di^m are the remarkable integers in the title. These rational coefficients 
di^m appeared in the evaluation of the quartic integral 



(1.2) No4a;m) := 
for a > —1, m E N. The formula 

(1.3) NoA{a;m) 

with 



dx 



(a;4 + 2ax^ + l)™+i ' 
TT Pm{a) 



[2(a + 1)] 



rn+7 



2 - ■ 1 

\ " , J 



(1.4) Pmia) = } di^ma' 



1=0 



has been established by a variety of methods, some of which are reviewed in [3]. 
The symbolic status of p.2[) has not changed much since we last reported on [26] . 
Mathematica 6.0 is unable to compute it when a and m are entered as parameters. 
On the other hand, the corresponding indefinite integral is evaluated in terms of 
the AppcU-Fl function defined by 

(1.5) Fi{a;bi,b2;c;x,y) := ^ ^ — x y 

m.n.[c)m+n 
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+ 2aa;2 + 1) 
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Here {a)k = a{a + 1) • • • (a + fc — 1) is the ascending 



where a± :— a ± V— 1 + 
factorial. 

The coefBcients {d/.m : < ^ < to} have remarkable properties that will be 
discussed here. Those properties have mainly been discovered by following the 
methodology of Experimental Mathematics, as presented in [UlIT^. Many of the 
properties presented here have been guessed using a symbolic language and subse- 
quently established by traditional methods. The reader will find in [8] a detailed 
introduction to the polynomial Pmio) in (|l-4p . 



2. A TRIPLE SUM EXPRESSION FOR d, 



Our first approach to the evaluation of p.3|) was a byproduct of a new proof of 
Wallis's formula. 



(2.1) J2,m '■— 



dx f 2to 



where to is a nonnegative integer. Wallis' formula has the equivalent form 
TT 2 2 4 4 2n 2n 



(2.2) 



2 1 3 3 5 2n - 1 2n + 1 



(2.4) G(a,6) = 



The reader will find in [5] a proof of the equivalence of these two formulations. 

We describe in [9] our first proof of (|2.ip . Section [3] shows that a simple exten- 
sion leads naturally to the concept of rational Landen transformations. These are 
transformations on the coefficients of a rational integrand that preserve the value 
of the integral. It is the rational analog of the well known transformation 

(2.3) a — , b 1-^ \fab 

that preserves the elliptic integral 

-"/^ dd 

\/ cos^ 6* + 6^ sin^ Q 

The reader will find in [T3] and [Mj details about these topics. 

The proof of Wallis' formula begins with the change of variables x = tan0. This 
converts J2,m to its trigonometric form 

(2.5) J.^„^^\^\o.--BdB = 

The usual elementary proof of (|2.5p presented in textbooks is to produce a recur- 
rence for J2.m- Writing cos^ = 1 — sin^ Q and using integration by parts yields 

_ 2to - 1 

(2.6) J2,m — — Jl,m-\- 

Am 

Now verify that the right side of (|2.5p satisfies the same recursion and that both 
sides give 7r/2 for to = 0. 
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A second elementary proof of Wallis's formula, also given in [9j, is done using a 
simple double-angle trick: 

J2,m = y cos'^"^ 9 dO^ J ( 1 d9. 

Now introduce the change of variables ip = 29, expand and simplify the result by 
observing that the odd powers of cosine integrate to zero. Hence (|2.5p is reduced 
to an inductive proof of the binomial recurrence 

Lm/2J 

(2.7) J2,™ = 2-™ U H^.. 



1=0 



Note that J2,m is uniquely determined by (|2.7p along with the initial value J2,o 
7r/2. Thus (|2.5p now follows from the identity 

L™/2J / \ /o-\ /o 

(2.8) fim):= 
since (12.81) can be written as 



2i J \ i J \ m 



L™/2J 

i=0 



2i 



where 



)2i+l 



The last step is to verify the identity (|2.8p . This can be done mechanically using 
the theory developed by Wilf and Zeilberger, which is explained in [571 [28]. The 
sum in (|2.8p is the example used in [28] (page 113) to illustrate their method. 

Note. The WZ-method is an algorithm in Computational Algebra that, among 
other things, will produce for a hypergeometric/holonomic sum, such as (13.71) . a 
recurrence like (|3.10p . The reader will find in [27 and information about this 
algorithm. 

The command 

ct{binomial(m, 2i) binomial{2i, i)2^^*, 1, i, to, N) 

produces 

(2.9) /(to + 1) = _t^/(TO), 

TO + 1 

a recursion satisfied by the sum. One completes the proof by verifying that 2^™ (^^™) 
satisfies the same recursion. Note that (12. 6p and (|2.9p are equivalent since J2,m and 
/(to) differ only by a factor of 7r/2'"+^. 

We have seen that Wallis's formula can be proven by an angle- doubling trick 
followed by a hypergeometric sum evaluation. Perhaps the most interesting appli- 
cation of the double-angle trick is in the theory of rational Landen transformations. 
See [24l for an overview. 
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Now we employ the same ideas in the evaluation of (II. 3p . The change of variables 
X = tan 6 yields 

. X r^V cos^e \"+' de 

Nn 4{a;m) — / ^ ^ x 

' Jo \sm^9 + 2asm^ecos^9 + cos'^dJ cos^ 

Observe first that the denominator of the trigonometric function in the integrand 
is a polynomial in m = 26. In detail, 

sin" + 2a sin^ cos^ + cos" = 2 [(1 + a) + (1 - a) cos^ u] . 

In terms of the double-angle u — 20, the original integral becomes 



A'o,4(a; m 



_ r,-(m+l) 



(1 + cosu)^ 



(1 + a) + (1 — a) cos^ u J 1 + cos u 

Next, expand the binomial (1 + cosu)'^'"+^ and check that 

/•TT 

(2.10) / [(l + a) + (l-a)cos2u]^'™^^ 'cos^'wdM = 



for j odd. The vanishing of half of the terms in the binomial expansion turns out to 
be a crucial property. The remaining integrals, those with j even, can be simplified 
by using the double-angle trick one more time. The result is 

No4a; ™) = £ £ [(3 + a) + [1 - a) cos w]"(™+") (1 + cos vY dv, 

where v = 2u and we have used the symmetry of cosine about v = tt to reduce 
the integrals form [0, 27r] to [0,7r]. The familiar change of variables z — tan(v/2) 
produces (jl.3p with the complicated formula 

' ' (— 1)'^^'^'' /2fc\ /2m + 1\ /m — s — A /s + A /fc — s — j 



di,m ^Y.Y. 



^ 23fe \k)\2s^2i)\ m-k l\ j I \ l-j 

j=0 s=Q k=s+l \/\ J/ \ / \ J / \ J 

Note. In spite of its complexity, obtaining this expression was the first step in the 
mathematical road described in this paper. It was precisely what Kauers and Paule 
[20j required to clarify some combinatorial properties of di^m- Some arithmetical 
properties can be read directly from it. For example, we can see that di^m is a 
rational number and that i^^di^rn € that is, its denominator is a power of 2 
bounded above by 3m. Improvements on this bound are outlined in Section [31 

3. A SINGLE SUM EXPRESSION FOR d/.m 

The idea of doubling the angle that proved productive in Section [5] can be ex- 
pressed in the realm of rational functions via the change of variables 

-^2-1 



(3.1) V = R2{x) - 2^ . 

The inverse has two branches 

(3.2) x = y±^v^ + l, 

where the plus sign is valid for x £ (0, -l-oo) and the other one on (— oo,0). The 
rational function R2 arises from the identity 

(3.3) cot26' = i?2(cot6'). 
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This change of variables gives the proof of the next theorem. 

Theorem 3.1. Let f be a rational Junction and assume that the integral of f over 
R is finite. Then 



(3.4) 



fix) dx 



fiy + V7+1) + f{y - ViFTi)\ dy + 

ydy 



fiy + vV + i) - fiy - Vy' + 1) 



Moreover, if f is an even rational function, the identity Jg.^p remains valid if one 
replaces each interval of integration by M+ . 



Theorem 3.2. For to e N, let 



(3.5) 

Define 



Qix) 



1 



iiy) 



Qiy + V7TT) + Qiy - V^F+i) 
y 



Then 
(3.6) 
where 
(3.7) 



Qiiy) - 



Qiy + Vy' + 1) - Qiy - ViF+i) 

Tmi2y) 
2™(1 + a + 2y2)m+i' 

' TO + A;^ 



Proof. Introduce the variable cj) — y + yy^+T. Then y — y^y'^ + 1 = —4> ^ and 
y = i(0 — (^~^). Moreover, 

Qiiy) = [Q(0) + g(r')] +|J^(Q(0)-Q(r')) 

:= S'm(0)- 

The result of the theorem is therefore equivalent to 

(3.8) 2™ (1 + a + i(0 - r ')'))'"^' ^™(0) = r™(<^ - r')- 

A direct simplification of the left hand side of (|3.8p shows that this identity is 
equivalent to proving 



(3.9) 



■1+1 + <^-(2m+l) 
0+0-1 



— T 

J- 77 



To establish this, one simply checks that both sides of 
order recurrence 



satisfy the second 



(3.10) 



Cm+2 



^)c„i+l + Cm = 0, 
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and the values for m ~ and m = 1 match. This is straight-forward for the 
expression on the left hand side, while the WZ-method settles the right hand side. 

□ 



We now prove (|1.3p . The identity in Theorem 13.11 shows that 

/•OO /"OO 

(3.11) / Q{x)dx^ / Qi{y)dy, 

Jo Jo 

and this last integral can be evaluated in elementary terms. Indeed, 



2y2) 



2\m+l 



^_^(m + k\ [°° {2yf^ dy 



2™ ^ Vm - fcy Jo (1 + a + 2y2)m+i ' 



The change of variables y ~t \/l + gives 

\m, — 

k=0 

and the elementary identity 



(l+f2)m+l^ 



i^*: TT /2fc\ /2m - 2fc\ /m^ ^ 



(1 + 22"+i 



gives 



fe=0 

This can be simplified further using 



m + k\ / 2k\ / m + k\ / m 
m — k I \ k J \ rn J \ k 



(3.12) 

and the equality (|3.1ip to produce 

(^•1^) L Q(^)-^^-2^ [2(i+«)r+v2 g/^( ^ )(m-0(i+«)'=- 

This completes the proof of (|1.3p . The coefficients d/.^ are given by 

m 

(3.14) = 2-2'"^2* 



k=l 



. /2to - 2A:\ /m + /c\ //c 
Vm — /c/V m J \ l 



This is clearly an improvement over the expression for di^m given in the previous 
section. 

We now see that di^m is a positive rational number. The bound on the denomi- 
nator is now improved to 2m — 1. This comes directly from p.l4p and the familiar 
fact that the central binomial coefficients f^"*) are even. 
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4. A FINITE SUM 

The previous two sections have provided two expressions for the polynomial 
Pm{a)- The elementary evaluation in Section [5] gives 

(4.1) 

and the results described in Section [3] provide the alternative expression 
(4.2) P„(a) = 2-'"f^2-'^-n^)f^™"'=)(a + ir-^ 



fe=0 



The reader will find details in [5]. Comparing the values at a = 1 given by both 
expressions leads to 

i:-"(!)C';r)^s-^'(?)C";,:' 

The identity (|4.3p can be verified using D. Zeilberger's package EKHAD [55]. In- 
deed, EKHAD tells us that both sides of (|4.3p satisfy the recursion 

(2m + 3) (2m + 2)/(m + 1) = (4m + 5) (4m + 3)/(m). 

To conclude the proof by recursion, we check that they agree at m = 1 . A symbolic 
evaluation of both sides of (|4.3p leads to 

^^^^ 22™+ir(2TO + 3/2) _ 22™+i0F 



0Fr(2m + 2) r(-2m- l/2)r(2TO- 

The identity (|4.3p now follows from 



(4.5) r(m+i) = ^(Mlfor^eN. 

An elementary proof of (|4.3|) would be desirable. 

The left hand sum admits a combinatorial interpretation: multiply by 2^™"*"^ to 
produce 

,«) s.C")-!: 

Consider the set X of all paths in the plane that start at (0, 0) and take 2m + 1 
steps in any of the four compass directions {N — (0,1), S = (0, —1), E = (1, 0) and 
W = (—1,0)) so that the path ends on the y-axis. Clearly there must be the same 
number of E's and W's, say j of them. Then to produce one of these paths, choose 
which is E and which is W in ( j*) ways. Finally, choose the remaining 2m + 1 — 2j 
steps to be either N or S, in 2^"*+^^^^ ways. This shows that the set X has 5*1 (m) 
elements. 

Now let Y be the set of all paths of the x-axis that start and end at 0, take steps 
e = 1 and w = — 1, and have length 4m + 2. The cardinality of Y is clearly (jj^+i) • 
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There is a simple bijection between the sets X and Y given hy E ^ ee, W 
WW, N —^ ew, S —^ we. Therefore, 



(4.7) Slim) 



/4to + 2 
\2m+ 1 



We have been unable to produce a combinatorial proof for the right hand side 
of (lOl). 

5. A RELATED FAMILY OF POLYNOMIALS 

The expression p.l4p provides an efficient formula for the evaluation of di ^ 
when I is close to m. For example, 



(5.1) drn,,n = T 



and = (2m + l)2-(-+i) P"" 

/ V "1 



Our attempt to produce a similar formula for small I led us into a surprising family 
of polynomials. 

The original idea is very simple: start with 

and compute di^m as coming from the Taylor expansion at a = of the right hand 
side. This yields 

(5.3) d,™ = J——, ai{m) \{{Ak - I) - Pi{m)\{{Ak + I)] , 



k=l k=l 



where ai and (3i are polynomial in m of degrees I and ^ — 1, respectively. The explicit 
expressions 

L'/2J / , N m+t m t-1 

(5.4) «,(m)= ^ (^J n (4^-1) n (2^+1) 11(4'^+ 1)' 
and 

L(i+1)/2J , , m+t-1 m t-1 

(5.5) A(m)= Lt-l) n (4^+1) n (2^+1) 11(4'^- 1)' 
are given in [TU] , 

Trying to obtain more information about a; and /3; directly from (j5.4[ 15. 5p proved 
difRcult. One uninspired day, we decided to compute their roots numerically. We 
were pleasantly surprised to discover the following property. 

Theorem 5.1. For all I > 1, all the roots of ai{m) = lie on the line Rem = — ^• 
Similarly, the roots of f3i {m) = for I > 2 lie on the same vertical line. 

The proof of this theorem, due to J. Little !23], starts by writing 

(5.6) Ms) aiiis - l)/2) and Bi{s) := A((s - l)/2) 
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and proving that Ai is equal to II times the coefficient of it' in /(s, u)g{s, u), where 
/(s, u) = {1 + 2^)"/^ and g{s, u) is the hypergeometric series 

. X . X fs 111 

(5.7) gis,u)^2Fii- + -,-;-;4u' 



2 4' 4' 2' 

A similar expression is obtained for Bi{s). From here it follows that Ai and Bi each 
satisfy the three-term recurrence 

(5.8) xi+i{s) = 2sxi{s) - (s2 - (21 - l)^)xi-i{s). 

Little then establishes a version of Sturm's theorem to prove the final result. 

The location of the zeros of a;(m) now suggest to study the behavior of this 
family as Z ^ oo. In the best of all worlds, one will obtain an analytic function of 
m with all the zeros on a vertical line. Perhaps some Number Theory will enter 
and ... one never knows. 

6. Arithmetical properties 
The expression (|5.3[) gives 

m m 

(6.1) m!2™+i di,™ = (2771 + 1) [| (4fc - 1) - (4fc + 1), 

k=l k=l 

from where it follows that the right hand side is an even number. This led naturally 
to the problem of determining the 2-adic valuation of 

ni m 

(6.2) A,,™ := ;!m!2"+'d,,™ = aiim)Y[{4k - 1) ~ Piim)Y[{4k + 1) 

Recall that, for a; G N, the 2-adic valuation 1^2(2^) is the highest power of 2 that 
divides x. This is extended to x = a/6 G Q via 1^2 (a^) ~ 1^2(0,) — 7^2(6), leaving 7^2(0) 
as undefined. It follows from (I6.3|l that 

(6.4) = 2"(2m)! and = 2"-i(2m - l)!(2m + 1), 

so these 2-adic valuations can be computed directly from Legendre's classical for- 
mula 

(6.5) 1^2(2;) = X - S2{x), 

where S2{x) counts the number of I's in the binary expansion of x. 
At the other end of the /-axis, 

m 

(6.6) Ao,m = l[{4k - 1) 

k=l 

is clearly odd, so i^2(^o,t?i) — 0. The first interesting case is / = 1: 

m rn 

(6.7) Ai,™ - (2m + l)l[{4k-l)-l[{4k + l). 

k=l fe=l 

The main result of [10] is that 

(6.8) V2{Al,rn) = V2im{m + 1)) + 1. 
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This was extended in [2]. 
Theorem 6.1. The 2-adic valuation of Ai^„i satisfies 

(6.9) J^2(A;,„) = Z/2((to + 1-020+^ 

where {a)k = a(a + 1) • • • (a + fe — 1) is the Pochhammer symbol for k > 1. For 
k — 0, we define (a)o — 1. 

The proof is an elementary application of the WZ-method. Define the numbers 

(6-10) S;.™ := 1''"^ n ' 

2'(m + 1 - 1)21 

and use the WZ-method to obtain the recurrence 

Bl-l^rn = (2m + l)Bi^rn - (m - l){m + I + l)Bi + i^,n, I <l <m - I. 

Since the initial values i?m,m — 1 and -Bm-i,m — 2m + 1 are odd, it follows induc- 
tively that m is an odd integer. The reader will also find in 2 a WZ-free proof 
of the theorem. 

Note. The reader will find in [3] a study of the 2-adic valuation of the Stirling 
numbers. This study was motivated by the results described in this section. The 
papers [HI [HI [22l [211 [33l 34] contain information about 2-adic valuations of related 
sequences. 



7. The combinatorics of the valuations 

The sequence of valuations {1^2(^1, m) m > 1} increase in complexity with I. 
Some of the combinatorial nature of this sequence is described next. The first 
feature of this sequence is that it has a block structure, reminiscent of the simple 
functions of Real Analysis. 

Definition 7.1. Let s € N, s > 2. We say that a sequence {aj : j G N} has block 
structure if there is an s G N such that each t G {0, 1, 2, • • • }, we have 

(7.1) ast+i = ast+2 — ■ ■ ■ — a5(t_|_i). 

The sequence is called s-simple if s is the largest value for which (|7.ip occurs. 

Theorem 7.2. For each I > 1, the set X{1) :— {v2{Ai.m) : m> 1} is an s-simple 
sequence, with s = 2^^'^^^''\ 

We now provide a combinatorial interpretation for X{1). This requires the maps 

F{{ai, a2, as, ■ ■ ■}) := {ai, ai, 02, 03, • • • }, 

T{{ai, a2, as, ■ ■ ■}) := {ai, 03, 05, 07, • • • }. 

We will also employ the notation c {z^2(m) : m > 1} = {0, 1, 0, 2, 0, 1, • • • }. 

We describe an algorithm that reduces the sequence X{1) to a constant sequence. 
The algorithm starts with the sequence X{1) := {i'2{Aij+rn-i) ■ m > 1 } and then 
finds and n G N so that X{1) is is 2"-simple. Define ¥(1) := T" {X{1)). At the 
initial stage. Theorem 17. 21 ensures that n = 1'2(0- The next step is to introduce 
the shift Z{1) := Y{1) - c and finally define W{1) := F{Z{1)). If W{1) is a constant 
sequence, then STOP; otherwise repeat the process with W instead of X. Define 
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Xk{l) as the new sequence at the end of the (fc — l)th cycle of this process, with 
Xi{l)^X{l). 

This algorithm produces a sequence of integers rij, so that Xk{l) is 2"*^ -simple. 
The integer vector := {rii, n2, n^, ■ ■ ■ is called the reduction sequence 

of I. The number ui is the number of cycles requires to obtain a constant sequence. 

Definition 7.3. Let I G N. The composition of I, denoted by is an integer 

sequence defined as follows: Write I in binary form. Read the digits from right 
to left. The first part of is the number of digits up to and including the 

first 1 read in the corresponding binary sequence; the second one is the number of 
additional digits up to and including the second 1 read, and so on until the number 
has been read completely. 

Tiieorem 7.4. Let {ki, ■ ■ ■ ,kn : < ki < k2 < ■ ■ ■ < kn}, be the unique collection 
of distinct nonnegative integers such that I ~ X^iLi ■ Then the reduction sequence 
fl{l) of I is {ki + 1, fc2 - fci, • • • , A;„ - fc„_i}. 

It follows that the reduction sequence n{l) is precisely the sequence of compo- 
sitions of I, that is, n(l) — ili{l). This is the combinatorial interpretation of the 
algorithm used to reduce X(l) to a constant sequence. 

8. Valuation patterns encoded in binary trees 

In this section we describe the precise structure of the graph of the sequence 
{i'2{A-i,m),T^^' > I}- The reader is referred to [3D] for complete details. In view 
of the block structure described in the previous section, it suffices to consider the 
sequences {i'2(C/,m), nT- > I}, which are defined by 

The emerging patterns are still very complicated. For instance, Figure [T] shows the 
case I = 13 and Figure [2] corresponds to I = 59. The remarkable fact is that in 
spite of the complexity of J^2(C'/,m) there is an exact formula for it. The rest of this 
section describes how to find it. 




Figure 1. The valuation 1/2(6*13,™) 



We describe now the decision tree associated to the index /. Start with a root vq 
at level fc = 0. To this vertex we attach the sequence {1^2(6';, m) : m > 1} and ask 
whether V2{Ci^m) — V2(m) has a constant value independent of m. If the answer is 




yes, we say that vq is a terminal vertex and label it with this constant. The tree is 
complete. If the answer is negative, we split the integers modulo 2 and produce two 
new vertices, wi, V2, connected to vq and attach to the classes {i^2(C'i.2m-i) : m > 1} 
and {i^2(C'i,2m) : TO > 1} to these vertices. We now ask whether !^2(C'i,2m-i)~'^2(TO) 
is independent of m and the same for 1^2 (C;, 2m) ~ J^2(™)- Each vertex that yields 
a positive answer is considered terminal and the corresponding constant value is 
attached to it. Every vertex with a negative answer produces two new ones at the 
next level. 

Assume the vertex v corresponding to the sequence {2^{m — 1) + a : m > 1} 
produces a negative answer. Then it splits in the next generation into two vertices 
corresponding to the sequences {2*^+^ (m—l)+a : to > 1} and {2''+^(to— l) + 2'''+a : 
TO > 1}. For instance, in Figure[3l the vertex corresponding to {4to : m > 1}, that 
is not terminal, splits into {8to : to > 1} and {8to — 4 : to > 1}. These two edges 
lead to terminal vertices. Theorem 18.11 shows that this process ends in a finite 
number of steps. 




Theorem 8.1. Let I G N and T{1) he its decision tree. Define k*{l) := [logj 
Then 
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1) T(l) depends only on the odd part of I; that is, for r S N, we have T{1) — TiT'l), 
up to the labels. 

2) The generations of the tree are labelled starting at 0; that is, the root is genera- 
tion 0. Then, for < k < k*{l), the k-th generation ofT{l) has 2^ vertices. Up to 
that point, T{1) is a complete binary tree. 

3) The k*-th generation contains 2^ — I terminal vertices. The constants asso- 
ciated with these vertices are given by the following algorithm. Define 

(8.1) jiil,k,a):=-l + 2{l + 2'' -a), 
and 

(8.2) j^(^lk,a)^l + k + l + iy2 {{ji + I - 1)!) + 1^2 {{I - Ji)!) • 
Then, for 1 < a < 2'^*+^ — I, we have 

(8.3) U2 (C,_2'=(m-l)+a) = z^2(m) + 7i (Z, fc, o). 

Thus, the vertices at the k*-th generation have constants given by ^\{l,k,a). 

4) The remaining terminal vertices of the tree T{1) appear in the next generation. 
There are 2(Z — 2*^*''^) of them. The constants attached to these vertices are defined 
as follows: let 

(8.4) j2{l, k, a) := -I + 2(1 + 2''+^ - a), 
and 

(8.5) h{l,k,a):=j2{l,k,a + 2''). 
Define 

(8.6) 72(/, k,a):=l + k-h2 + iy2 {{32 + I - 1)!) + V2 {{I - J2)!) , 
and 

(8.7) 73(/, k,a):=l + k + 2 + V2 {{h + I - 1)0 + 1^2 {{I - is)!) . 
Then, for 2'=*(')+i - Z + 1 < a < 2'=*('), we have 

(8.8) U2 (C;,2fe*(0+i(™-i)+a) = z^2(m) +72(Z,fc*(Z),a), 

and 

(8.9) 1^2 (Q,2'=*(i) + i(m-l)+Q+2'=*(i)) = i^2(TO) + jsH , k* (l) , a) , 

give the constants attached to these remaining terminal vertices. 

Wc now use the theorem to produce a formula for J^2(C'3,m). The value fc*(3) = 1 
shows that the first level contains 2^+^ — 3 = 1 terminal vertex. This corresponds 
to the sequence 2m — 1 and has constant value 7, thus, 

(8.10) U2 (C3,2m-l) = 7. 

The next level has 2(3 — 2^) = 2 terminal vertices. These correspond to the se- 
quences 4m and 4to— 2, with constant values 9 for both of them. This tree produces 

r7 + i^2(^) ifm=lmod2, 

(8.11) ,.2(^3,™) = i 9 + 1/2 (f) ifm = 0mod4, 

[9 + 1/2 (^) ifm = 2mod4. 



14 



VICTOR H. MOLL AND DANTE V. MANNA 



The complexity of the graph for I 
this valuation. The theorem yields 



13 is reflected in the analytic formula for 



.12) 



1^2 (Cl3,m) = < 



36 + V2 


(^) 


if 


m 




1 mod 8, 


37 + V2 


I 8 ; 


if 


m 




2 mod 8, 


36 + V2 


\ 8 ) 


if 
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3 mod 8, 


40 + vi 


\ 16 y 


if 


m 




4 mod 16, 


38 + J/2 


\ 16 / 


if 


m 




5 mod 16, 


39 + V2 


(^) 


if 


m 




6 mod 16, 


38 + Z/2 




if 


m 




7 mod 16, 


40 + vi 




if 


m 




8 mod 16, 


40 + vi 




if 


m 




12 mod 16, 


38 + V2 




if 


m 




13 mod 16, 


39 + V2 




if 


m 




14 mod 16, 


38 + ^2 


{^) 


if 


m 




15 mod 16, 


^40 + 1^2 




if 


m 




16 mod 16. 



The details for Theorem 18.11 are given in 



Note. The p-adic valuations of A;.„j for -p odd present phenomena different from 
those explained for the case p = 2. Figure 2] shows the plot of 1^17(^1, ,„) where we 
observe linear growth. Experimental data suggest that, for any odd prime p, one 
has 

vn 

(8.13) V-p[Ai^r, ^ 



Figure [5] depicts the error term j^i7(v4i, 
remains to be explored. 



p-1 
) — m/16. 



The structure of the error 




Figure 4. The valuation vx-j{Ax^ra) 



9. UnIMODALITY AND LOG-CONCAVITY 



A finite sequence of real numbers {oq, ai, • • • 
exists an index < j < m such that ao < ai < 



1} is said to be unimodal if there 



< Gj and aj 



> a 



j+i 



> 



> an 
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A polynomial is said to be unimodal if its sequence of coefficients is unimodal. The 
sequence {ao, ai, • • • , a„i} with Oj > is said to be logarithmically concave (or log- 
concave for short) if aj+iaj^i < a'j for 1 < J < m — 1. It is easy to see that if a 
sequence is log-concave then it is unimodal [36j . 

Unimodal polynomials arise often in combinatorics, geometry, and algebra, and 
have been the subject of considerable research in recent years. The reader is referred 
to |29j and [14j for surveys of the diverse techniques employed to prove that specific 
families of polynomials are unimodal. 

For m € N, the sequence {di^m ■ < I < m} is unimodal. This is a consequence 
of the following criterion established in [6]. 

Theorem 9.1. Let ak be a nondecreasing sequence of positive numbers and let 
^i^) = X^fcLo o-kX^ ■ Then A{x + 1) is unimodal. 

We applied this theorem to the polynomial 
(9.1) j( ^ 

that satisfies Pm{x) = A{x + 1). The criterion was extended in to include the 
shifts A(x+j) and in [32] for arbitrary shifts. The original proof of the unimodality 
of Pm{o) can be found in |7]. 

In [26] we conjectured the log-concavity of {di^„i : < ^ < to}. This turned out 
a more difficult question. Here we describe some of our failed attempts. 

1) A result of Brenti [14] states that if A{x) is log-concave then so is A{x + 1). 
Unfortunately this does not apply in our case since (|9.ip is not log-concave. Indeed, 

2 

X 



-)4m-2*: (2 



2m \ /m + k 
m — k J \ TO 



^ k{m - k){2m - 2k + l){m + k + 1) 



[k + l)(m + fc)(2m - 2fc - l)(m - fc + 1)^ 

and this last factor could be negative — for example, for m — b and j — 4. The num- 
ber of negative terms in this sequence is small, so perhaps there is a way out of this. 
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2) The coefficients di^m satisfy many recurrences. For example, 

, s 2m + 1 , , , (m + j){ni + 1 — ?') , , , 

(9.2) d,+i(m) = —-^d,{m) - ^ ^^i^^^ ^d,_i(m). 

This can be found by a direct appUcation of WZ method. Therefore, di^m is log- 
concave provided 

(9.3) ^{2^^ + l)dj^i{m)dj{m) < {m + j){m + 1 - j)dj^i{m)'^ + j{j + l)dj{m)'^ . 
We have conjectured that the smallest value of the expression 

(9.4) (m + j)(m + 1 - j)dj^i{mf + j{j + l)dj(m)^ - j{2m + l)dj^i{m)dj{m) 

is 2^™TO(m + and it occurs at j = to. This would imply the log-concavity 

of {d/,m '■ < I < m}. Unfortunately, it has not yet been proven. 

Actually we have conjectured that the d/^m satisfy a stronger version of log- 
concavity. Given a sequence {aj} of positive numbers, define a map 

by bj := Oj — aj^iaj^i. Thus {aj} is log-concave if {bj} has positive coefficients. 
The nonnegative sequence {aj} is called infinitely log- concave if any number of ap- 
plications of £ produces a nonnegative sequence. 



Conjecture 9.2. For each fixed m G N, the sequence {di^m : < ^ < to} is 
infinitely log-concave. 

The log-concavity of {di^m : < ^ < has recently been established by M. 
Kauers and P. Paule '^Hj as an applications of their work on establishing inequalities 
by automatic means. The starting point is the triple sum expression in Section [2] 
written as 

_ (-1)'=+J-Y2to+1\ /to-s\ /2(fc + s)\ /s\ / k 
h^^^^ 2. )[ k )[ k + s )[j)[l-j 

Using the RISC package Multisum [3S] they derive the recurrence 

(9.5) 2(m -I- l)d;,™+i = 2(1 + m)di_i,™ -I- {21 + 4m + 3)^;,™. 

The positivity of di^m follows directly from here. To establish the log-concavity of 
di^rn the new recurrence 

Al{l + l)di+i^,n = -2(2/ - 4to - 3)(/ + m + l)rf/,„i + 4(/ - to - 1)(to + l)di,„i+i 

is derived automatically and the log-concavity of di^m is reduced to establishing the 
inequality 

^ 4(m + 1) (4((/ - m - 1)(to + 1) - {2P - 4to^ -7m- 2.)di^m+idi,m) 
- 16to3 + 16/m2 + 40m^ -f 28/m -f- 33to + 91 + 9 

The 2-log-concavity of {d;,™ : < / < m}, that is {{di^m}) > {0,0, ...0} re- 
mains an open question. At the end of [20] the authors state that "...we have little 
hope that a proof of 2-logconcavity could be completed along these lines, not to 
mention that a human reader would have a hard time digesting it." 
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The general concept of infinite log-concavity has generated some interest. D. 
Uminsky and K. Yeats [31] have studied the action of £ on sequences of the form 

(9.6) {••• ,0,0, l,xo,xi, • • • ,a;„,--- , xi, a;o, 1, 0, 0, • • • } 

and 



(9.7) 



{••• ,0,0, l,a;o,a;i,--- ,a;„,a;„,--- , xi, xq, 1, 0, 0, • • • } 



and established the existence of a large unbounded region in the positive orthant 
of R" that consists only of infinitely log-concave sequences {xq, . . . ,Xn}- P. McNa- 
mara and B. Sagan [25] have considered sequences satisfying the condition a| > 
rofc_iafc_|_i. Clearly this implies log-concavity of r > 1. Their techniques apply to 
the rows of the Pascal triangle. Choosing appropriate r-factors and a computer 
verification procedure, they obtain the following. 

Theorem 9.3. For fixed n < 1450, the sequence {(^') : < k < n} is infinite 
log-concave. 

In particular, they looked for values of r for which the r-factor condition is preserved 
by the £-operator. The factor that works is r = ^'^./^ (the square of the golden 
mean). This technique can be used on a variety of finite sequences. See |25j for a 
complete discussion of the techinque. 

McNamara and Sagan have also considered (/-analogues of the binomial coeffi- 
cients. In order to describe these extensions we introduce the basic notation and 
refer to the reader to [19] and [5] for more details on the world of ^-analogues. 

Let g be a variable and for n G N define 

(9.8) [n], = ^-^ = l + q + q^ + --- + q--\ 

The Gaussian-polynomial or q-binomial coefficients are defined by 

\r,]\ 



(9.9) 



The Gaussian polynomials have nonnegative coef- 



where [n]q\ := [l],^], • • • [i 
ficients. We will say that the sequence of polynomials {fkiq)} is q-log-concave if 
£(/fc(g)) is a sequence of polynomials with nonnegative coefficients. The extension 
of this definition to infinite q-log- concavity is made in the obvious way. 
Observe that 

(9.10) 

McNamara and Sagan have established the surprising result: 



is not infinite q-logconcave. 



In fact they established that applying £ twice gives polyomials with some nega- 
tive coefficients. As a compensation, they propose: 





— lim 


n 






9^1 


k 







n 




sequence < 




k 


i 



Conjecture 9.5. The sequence 
k>0. 




is infinite q-log-concave for all fixed 



1 ) n>k 
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Another g-analog of the binomial coefficients that arisies in the study of quantum 
groups is defined by 

(9.11) (n) := ^— = ^(1 + q^ + q^ + ... + q^-\ 

From here we proceed as in the case of Gaussian polynomials and define 

/ n\ {n)\ 



^^'^'^^ {k)\{n-k)\ 

where {n)\ = (1)(2) ■ ■ ■ (n). For these coefficients McNamara and Sagan have pro- 
posed 

Conjecture 9.6. a) The row sequence {(fe)}j,>Q infinitely q-log-concave for all 
n>0. 

b) The column sequence {{2}}n>k infinitely q-log-concave for all fixed n>0. 

c) For all integers < u < v, the sequence {("m™")}„j>o infinitely q-log-concave 
for all n > 0. 

This conjecture has been verified for all n < 24 with v < 10. When u > v, using 



(9.13) 



k I (/"'^ 



one checks that the lowest degree of ("^")^ — ("^u") ^'^ ^'^'^ sequence is not 

even g- log-concave. Sagan and McNamara observe that when u = v, the quantum 
groups analoge has exactly the same behavior as the Gaussian polynomials. 

Newton began the study of log-concave sequences by establishing the following 
result (paraphrased in Section 2.2 of [18]). 

Theorem 9.7. Let {ak\ he a finite sequence of positive real numbers. Assume all 
the roots of the polynomial 

(9.14) P[ak] x] := ao + aix H h a„a;" 

are real. Then the sequence {ak\ is log-concave. 

McNamara and Sagan [^S] and, independently, R. Stanley have proposed the 
next conjecture. 

Conjecture 9.8. Let {a^} be a finite sequence of positive real numbers. If P[ak]x\ 
has only real roots then the same is true for P[£(afc);x]. 

This conjecture was also independently made by Fisk. See [25j for the complete 
details on the conjecture. 

The polynomials Pm{o) in (|1.4p are the generating function for the sequence 
{di^m} described here. It is an unfortunate fact that they do not have real roots [7] 
so these conjecture would not imply Coniecture l9.2l In spite of this, the asymptotic 
behavior of these zeros has remarkable properties. Dimitrov [17] has shown that, 
in the right scale, the zeros converge to a lemniscate. 

The infinite-log-concavity of {dj.m} has resisted all our efforts. It remains to be 
established. 
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